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ABSTRACT 

We present evolutionary models for young low mass stars and brown dwarfs taking into 
account episodic phases of accretion at early stages of the evolution, a scenario supported by 
recent large surveys of embedded protostars. An evolution including short episodes of vigorous 
accretion (M > lCP 4 M yr -1 ) followed by longer quiescent phases (M < 10 _6 MQyr _1 ) can 
explain the observed luminosity spread in HR diagrams of star forming regions at ages of a 
few Myr, for objects ranging from a few Jupiter masses to a few tenths of a solar mass. The 
gravitational contraction of these accreting objects strongly departs from the standard Hayashi 
track at constant T e fj. The best agreement with the observed luminosity scatter is obtained if 
most of the accretion shock energy is radiated away. The obtained luminosity spread at 1 Myr 
in the HR diagram is equivalent to what can be misinterpreted as a ~ 10 Myr age spread for 
non-accreting objects. We also predict a significant spread in radius at a given T e ff, as suggested 
by recent observations. These calculations bear important consequences on our understanding of 
star formation and early stages of evolution and on the determination of the IMF for young (< 
a few Myr) clusters. Our results also show that the concept of a stellar birthlinc for low-mass 
objects has no valid support. 

Subject headings: stars: formation — stars: low-mass, brown dwarfs — accretion, accretion disks 



1. Introduction 

The significant luminosity spread observed in 
Hertzsprung-Russell (HR) diagrams of star form- 
ing regions (SFRs) and young clusters is a well 
known feature, which has been confirmed with the 
improvement of observational techniques (see Hil- 
lenbrand 2008 and references therein). Whether 
this luminosity spread arises from a physical pro- 
cess, observational uncertainties or reveals a sig- 
nificant age spread is a crucial question, with im- 



portant consequences on our understanding of star 
formation (Hartmann 2001). This hypothetical 
age spread is used as an argument in favour of 
a slow star formation process, in conflict with 
other observational constraints and our current 
understanding of star formation (Hartman 2001; 
Ballesteros-Paredes & Hartmann 2007). Moti- 
vated by this controversy, we have conducted a 
systematic analysis of how accretion affects the 
evolution of young low mass stars (LMS) and 
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brown dwarfs (BDs) in order to explore the sensi- 
tivity of evolutionary tracks to the early accretion 
history. In this Letter, we present the first consis- 
tent evolutionary models for young LMS and BDs 
taking into account non-steady accretion phases 
at very early stages of the evolution. We show 
that this scenario can explain at least partly the 
observed luminosity spread in HR diagrams, with- 
out invoking an age spread. 

Current observational analysis of embedded 
protostars strongly suggest that accretion onto 
forming stars must be transient, with very large 
fluctuations (Dunham et al. 2008; Enoch et al. 
2009; Evans et al. 2009). Enoch et al. (2009) find 
in three clouds a large population of low luminos- 
ity class I sources that aggravate the well known 
"luminosity problem" (Kenyon et al. 1990). Long 
quiescent phases of accretion (M <, lO _6 M0yr _1 ) 
interrupted by short episodes of high accretion 
(M > lO _5 M0yr _1 ) provide a consistent picture 
explaining both the large population of low lumi- 
nosity class I sources and the small fraction of very 
luminous sources. Enoch et al (2009) also rule out 
drastic changes in the accretion rates from class 

to class I, and in particular the standard pic- 
ture of short (~ 10 4 yr) class duratiorQ Evans 
et al. (2009) suggest that a star could assemble 
half its mass during a few episodes of high accre- 
tion, occurring throughout about 7% of the class 

1 lifetime. The idea of non steady, time vary- 
ing accretion rates is not new, since for decades 
FU Ori objects have been providing evidence for 
short episodes of rapid accretion at early stages of 
evolution, with rates much larger than the afore- 
mentioned typical class-0 infall rates for low mass 
objects (Kenyon & Hartmann 1995). 

The calculations presented in this Letter are 
a theoretical formulation of such an evolution in- 
cluding phases of episodic accretion for proto- or 
young LMS and BDs. In §2, we briefly summa- 
rize the evolutionary models and the treatment of 
accretion; details will be presented in a forthcom- 
ing paper (Gallardo et al. 2009). Results and 
comparison with observations are presented in §3, 
followed by discussion and conclusion in §4. 



1 Standard estimates for ages and rates are ~ 10 4 yr and 
M ~ 10~ 5 MQyr- 1 for class and ~ 10 5 yr and M ~ 
10 — 6 M0yr — 1 for class I sources. 



2. Evolutionary models and treatment of 
accretion 

The evolutionary calculations for LMS and BDs 
are based on the Lyon stellar evolution code with 
input physics described in Chabrier & Baraffc 
(1997) and Baraffe et al. (1998). The treatment 
of accretion is based on a simplified ID approach, 
adopting similar assumptions and simplifications 
as Hartmann et al. (1997) and Siess et al. (1997). 
We assume that accretion onto the central object 
rapidly proceeds non-spherically, affecting only a 
small fraction 8 of the contracting object's sur- 
face. The object can thus freely radiate its energy 
over most of its photosphere (see Hartmann et al. 
1997 and references therein). The accreting mate- 
rial brings, per unit mass, a gravitational energy 
—GM/R and an internal energy +eGM/R, i.e. an 
energy rate: 



dt 



(e-1) 



GMM 
R 



(1) 



The value of e depends on the details of the accre- 
tion process, with e < 1 for gravitationally bound 
material and e < \ if gas accretes from a thin disk 
at the object's equator (Hartmann et al. 1997). 
We denote a (a < 1) the fraction of accreting in- 
ternal energy absorbed by the proto-star/brown 
dwarf, which thus contributes to its heat content. 
The total additional energy rate gained by the ac- 
creting object and the accreting luminosity radi- 
ated away thus read, respectively, for 6 <C 1 (see 
e.g. Hartmann et al. 1997): 



GMM 

L a dd = ae „ ! L a 

K 



. GMM , . 
e(l-a)— — (2) 



The case a <C 1 corresponds to accreted matter 
arriving on the object's surface with a lower spe- 
cific entropy than the object's one. 

3. Results 

3.1. Exploring accretion rates 

We have conducted evolutionary calculations 
taking into account the effect of accretion on proto 
low mass objects for a wide range of initial proto- 
star/BD masses from 10 _3 M Q (~ 1 Afj up ) to 
0.1 Mq, with arbitrary large initial radii between 
~ li?0 and ~ 4i?0. The results are compared to 
recent surveys of LMS and BDs between ~ 0.01 
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Table 1: Parameters of the evolutionary sequences shown in Fig. [T] m, and mf are respectively the initial (in 
Mj up ) and final (in Mq) masses, T c m and L\ the effective temperature and log of the luminosity (in units 
of L Q ) respectively at 1 Myr, T C ^ CAH9S and ^BCAHgs thc vames predicted by the BCAH98 (non accreting) 
models for the mass m,f at 1 Myr, M the accretion rate (in M Q yr _1 ) applied during a time At (in yr) and 
a the fraction of accretion energy absorbed by the object (assuming e=l/2, see text). 



case 


m i 


mf 


T c m 


Lx 


yBCAH98 


£BCAH98 


M 


At 


Q 


A 


1 


0.05 


2860 


-2.15 


2844 


-1.73 


10" 5 


5 10 3 





B 


1 


0.1 


3157 


-1.63 


3001 


-1.16 


10~ 5 


10 4 





C 


1 


0.2 


3320 


-1.27 


3193 


-0.69 


10~ 5 


2 10 4 





D 


1 


0.5 


3756 


-0.75 


3426 


-0.28 


5 10~ 5 


10 4 





E 


1 


0.1 


3000 


-1.16 


3001 


-1.16 


10~ 5 


10 4 


0.2 


F 


1 


0.1 


2998 


-1.16 


3001 


-1.16 


10~ 5 


10 4 


1 


G 


5 


0.1 


3051 


-1.53 


3001 


-1.16 


10~ 5 


10 4 





H 


10 


0.21 


3221 


-1.20 


3201 


-0.67 


10~ 5 


2 10 4 





I 


50 


0.55 


3571 


-0.59 


3468 


-0.22 


10~ 5 


5 10 4 





J 


50 


1.05 


4102 


-0.32 


3840 


0.17 


5 10~ 5 


2 10 4 





K 


100 


1.1 


3961 


-0.03 


3870 


0.20 


5 10~ 5 


2 10 4 





L 


100 


1.85 


4677 


0.31 


4377 


0.60 


5 10~ 5 


3.5 10 4 






Mq and 1 Mq in various SFRs, with characteristic 
ages of a few Myr, namely Taurus and Chamaeleon 
I (Muzerolle et al. 2005), p Ophiucus (Gatti et 
al. 2006), Orion molecular cloud (Peterson et al. 
2008), cr-Orionis (Gatti et al. 2008). Given the 
loose constraints on accretion rates at early times, 
during the class and class I embedded phases, 
we have explored a wide range of mass accretion 
rates and time dependences, and considered differ- 
ent values of a. We have considered (1) constant 
rates, (2) exponentially time-decreasing rates and 
(3) rates obeying the empirical mass dependence 
M oc M 2 observed in young clusters (Mohanty 
et al. 2005, Herczeg & Hillenbrand 2008). The 
strong constraint of our calculations is to be con- 
sistent with the observed rate determinations at 
an age of a few Myr. 

We find that assuming initial accretion rates 
~ lO _6 M0yr _1 , as traditionally used in proto- 
stellar evolutionary models (Myers et al. 1998; 
Young & Evans 2005), produces too small a lu- 
minosity scatter in HR diagrams after a few Myr 
to explain the observed spread. More severe ef- 
fects are obtained when adopting higher accretion 
rates, M ~ (1 — 5) x 10~ 5 M Q yr _1 , during the 
first few 10 3 to 10 4 years of evolution, depending 
on the mass. Fig. [I] shows the effect of such high 
early accretion rates on the evolution of low-mass 
objects for different initial masses rrii. Evolution 



proceeds as time increases from the right to the 
left part of the HR diagram for a given track. The 
tracks of the accreting objects are displayed up to 
an age of 10 Myr with the locations at 1 Myr in- 
dicated by the solid squares. The evolutionary se- 
quences start from large initial radii and thus with 
short thermal timescales, tkh ~ G r L ~ 10 3 yr, of 
the order of or less than the accretion timescale 
r acc ~ M/M. Consequently, variations of the ini- 
tial radius by a factor 2-3 have no significant effect 
on the overall evolution, barely changing the final 
location at 1 Myr. The parameters of the various 
sequences are given in Table [T| 

For a=0, i.e. if all the accreting energy is ra- 
diated away, the aforementioned high accretion 
rates yield significantly smaller radii than the ones 
of non accreting objects of same mass and age. 
This stems from the fact that, as mass builds 
up, the object's thermal timescale tkh rapidly in- 
creases, and becomes much longer than the accret- 
ing timescale r acc , so that rapid increase of gravi- 
tational energy is the only possibility to adjust to 
the accreting energy flow. For example, the case B 
sequence (solid blue line in Figs, [iji and|2| starts 
its evolution with m., = 1 Mj up and tkh ~ 10 3 yr. 
At 10 yr, when accretion is arbitrarily stopped, 
the object has a mass 0.1 M Q , a radius ~ 0.5 Rq, 
a luminosity L ~ 2.5 x 10~ 2 L Q and a thermal 
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Fig. 1. — Effect of accretion on the evolution in a 
HRD for objects with different initial masses m\. 
The two long dashed (black) curves are the 1 Myr 
and 10 Myr isochrones of Baraffe et al. (1998) for 
non accreting models. The solid (black) circles on 
the rightmost Y-axis give the luminosity of non- 
accreting models at 1 Myr for the indicated masses 
(in Mq). The coloured curves show calculations 
for various accretion rates described in Table [I] 
The solid squares on each curve gives the position 
at 1 Myr. (a): Initial mass m; = lMj up and a=0. 
(b): initial mass m; = lAfj up and a / 0. Note 
that both sequences reach the same location at 
1 Myr. (c): Initial masses m\ > lMj up . Other 
symbols are observations in star forming regions 
(V Gatti et al. 2006; . Gatti et al. 2008; + 
Peterson et al. 2008; A Muzerolle et al. 2005). 



timescale tkh ~ 2 x 10 7 yr. Similar results are 
obtained for all the other sequences starting from 
larger initial masses. The radius of these accret- 
ing objects is thus already smaller after a few 10 4 
yr than the one of the non-accreting counterparts 
of same mass at an age of 1 Myr (see Fig. [2j. 
The objects will eventually slowly contract toward 
their location on the HR diagram at 1 Myr after 
the high accretion phase has been stopped, look- 
ing much fainter than the non-accreting 1 Myr old 
objects (see Fig. [T]). The spread in luminosity ob- 
tained in Fig. [ljthus essentially reflects a spread in 
radius, as illustrated in Fig. [2j Note that we find 
similar results if we apply smaller, non-zero ac- 
cretion rates, M < 10 _6 M Q yr _1 , after the strong 
accretion phase, while fulfilling the condition to 
recover typical observed rates at ~ 1 Myr. 

For the case a/0, i.e if some fraction of the 
accretion energy contributes to the proto-star/BD 
thermal content, the spread at 1 Myr is smaller 
than for the a=0 case. Assuming e—1/2 in Eq. 
(2) , characteristic of accretion from a thin disk (see 
§2), we find that for a <J 0.2, the contraction of the 
structure due to mass accretion is partly compen- 
sated by this energy input, yielding a less drastic 
contraction. The objects are still slightly fainter 
at ~ 1 Myr than their non-accreting counterparts, 
but the effect is not as significant as in the a=0 
case. For a >0.2, the effect of the extra accreting 
energy contribution to the protostar heat content 
becomes dominant and the evolution now mostly 
proceeds at higher luminosity (see Fig. [lj>) and 
larger radius than for non-accreting objects. This 
in turn implies shorter thermal timcscalcs, so that 
once high accretion rates are stopped or signifi- 
cantly decreased to recover the observed values at 
~ 1 Myr, the object quickly contracts and reaches 
a position at 1 Myr in the HR diagram very close 
to the non accreting location. In order to main- 
tain a high luminosity and radius after ~ 1 Myr, 
i.e a, location in the HR diagram above the 1 Myr 
non accreting isochrone, too large accretion rates 
(and large values of a) have to be maintained, in- 
consistent with observations of class II objects at 
this age. 

3.2. Episodic accretion 

The exploratory analysis presented in the pre- 
vious section demonstrates that the early accre- 
tion history still significantly affects the evolution 
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Fig. 2. — Evolution of the radius versus T c g for 
accreting objects of different initial masses. The 
colored curves correspond to different accreting se- 
quences described in Table [l] The other curves 
and symbols are the same as in Fig. [T] The solid 
(black) c ircles on the rightmost Y-axis give the 
radius of non-accreting models at 1 Myr for the 
indicated masses (in Mq). 



of contracting objects after a few Myr. In this 
section, we show that the effects of accretion are 
similar, whether the object accretes continuously 
an amount of mass AM at a constant rate M ~ 
10~ 5 M Q yr~ 1 during At ~ 10 4 yr, as examined in 
§3.1, or whether it accretes the same amount of 
mass during a succession of short episodes of high 
accretion rates, M > lO _4 M0yT -1 , interrupted 
by longer quiescent phases. As mentioned in |T| 
episodic accretion seems to provide the most con- 
sistent explanation for current observations of pro- 
tostars, while a short (~ 10 4 yr) accretion phase 
based on previous estimates of class lifetimes, 
with accretion rates significantly decreasing dur- 
ing the subsequent class I phase, seems to be ruled 



out (Enoch et al. 2009; Evans et al. 2009). Sev- 
eral theoretical scenarios have been suggested to 
explain phases of non-steady accretion (Kenyon & 
Hartmann 1995; Vorobyov & Basu 2005; Zhu et 
al. 2008; Tassis & Muschovias 2005). We adopt 
here a burst mode as suggested by Vorobyov & 
Basu (2005) to explore the effect of episodic ac- 
cretion on the structure of proto-stars/BDs. We 
assume an arbitrary number of bursts, A^urstj with 
typical accretion rate Mburst > 10~ 4 M Q yr~ 1 and 
duration Aiburst ~ 100 yr, interrupted by qui- 
escent phases of duration At qu i ct — 1000-5000 yr 
(see Vorobyov & Basu 2005). During the quies- 
cent phases, we adopt M—0 as a simplificatiorj^] 
We have explored two possibilities for the begin- 
ning of the burst phase (cases a and b described 
in Table |2|. The results are displayed in Fig. 
[3] The burst phases are characterised by evo- 
lutionary tracks which follow an erratic behavior 
with abrupt variations of L and T c g. Our calcu- 
lations show that, depending on ra- u N^ mst and 
A^burst, it is possible to populate the region in 
the HR diagram after ~ 1 Myr or less between 
the (non-accreting) 1 Myr and 10 Myr isochrones, 
producing a natural spread in luminosity. For val- 
ues of Mburst > 10~ 4 M Q yr _1 , the evolution is 
severely affected while for smaller rates, the struc- 
ture is only moderately affected and the object has 
time to relax once episodic accretion stops, having 
properties at ~ 1 Myr similar to the non accreting 
counterpart of same mass and age. This is illus- 
trated by case 4a (see Table |2| in Fig. [3] (cyan dot- 
ted curve), with Mburst = 10 _5 M Q yr _1 . Similar 
evolutionary properties are obtained adopting case 
a or 6 for the beginning of the burst phase. The 
duration of the quiescent phase Ai qu j e t is found to 
be inconsequential and can be increased from 10 3 
yr to ~ 10 4 yr without significant effects. This 
reflects the fact that the thermal timescale of the 
accreting object rapidly exceeds values 3> 10 4 yr. 
The contracting object has thus no time to relax 
to a larger radius state for its new mass if qui- 
escent phases last less than 10 4 yr. For typically 
10 to 100 burst episodes, this means that episodic 
accretion can last a few 10 5 yr, in agreement with 
recent revised estimates of class and class I life- 
times (Enoch et al. 2009; Evans et al. 2009). 



2 Adopting rates M < 10 6 provides the same qualitative 
results and does not affect our conclusions. 
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Table 2: Parameters of the evolutionary sequences assuming episodic accretion and shown in Fig. [3j raj, 

TOf, Teffl, L 1 , TJCAH98 &nd £BCAH98 haye the 

same meaning as in Table 1. A/bm-st is the accretion rate 
during the bursts (in M yr _1 ), Ai quict the quiescent phase duration (in units of 10 3 yr) and -/V burst the total 
number of bursts. In all calculations, the burst duration varies between ~ 75 and ~ 100 years. The label a 
or b indicates whether the bursts start at the beginning of evolution (case a) or after a first phase of constant 
accretion (M = lO _5 M0yr _1 during At = 10 4 yr), as suggested by Vorobyov & Basu (2005) (case 6). 



case m; m f T cS1 L x T c B ff CAH98 ^bcah98 Mw At N ^ 



la 


1 


0.1 


3151 


-2.00 


3001 


-1.16 


10" 4 


1 


10 


lb 


1 


0.2 


3373 


-1.25 


3193 


-0.69 


io~ 4 


1 


10 


21) 


1 


0.41 


3650 


-0.90 


3384 


-0.46 


io~ 4 


1 


30 


3b 


10 


0.38 


3466 


-1.05 


3335 


-0.43 


10~ 4 


1 


30 


4a 


10 


0.05 


2833 


-1.77 


2844 


-1.73 


10~ 5 


5 


50 


5a 


50 


1.35 


4379 


-0.16 


4152 


0.37 


io~ 4 


1 


150 


6a 


100 


1.16 


4039 


0.00 


3924 


0.23 


3 10~ 4 


1 


40 


7a 


100 


1.67 


4482 


0.19 


4317 


0.53 


3 10~ 4 


1 


60 



4. Discussion and conclusion 

The main results presented in §3 can be summa- 
rized as follows: (i) for a — 0, an accreting object 
has a more compact structure at ~ a few Myr, 
i.ea smaller radius and thus a smaller luminosity 
- looking older -, than the non-accreting counter- 
part of same mass and age; (ii) if a fraction of the 
accretion energy is absorbed by the protostar inte- 
rior (a ^ 0), this extra energy source partly com- 
pensates or, for a J>0.2 (assuming e=l/2), even 
dominates the contraction due to mass accretion, 
leading in the latter case to a larger radius than 
the non-accreting object of same mass and age. 
Our results show that a scenario based on early 
stages of non-steady accretion, characterized by 
short duration (a few 100 yr) episodes of vigor- 
ous accretion (M ~ 10 _4 M Q yr _1 ), followed by 
longer quiescent (M < 10 _5 MQyr _1 ) phases, nat- 
urally produces a spread in the HR diagram at 
ages of ~ a few Myr, providing non spherical ac- 
cretion (through disk or funnels) occurs on a va- 
riety of proto-star/BD initial masses, from a few 
Mj U p to a few tenths of Mq. This scenario can 
easily produce a luminosity spread equivalent to 
an age spread of ~ 10 Myr (for non-accreting ob- 
jects), even though the objects are just 1 or so Myr 
old. We find that such a spread can be obtained 
if only a negligible amount of the accretion energy 
contributes to the contracting object's internal en- 
ergy (a <C 1), i. e if most of the accreting kinetic 
energy is radiated away. 



Our scenario, however, cannot explain the ob- 
served population of very luminous low mass ob- 
jects, which lie largely above the 1 Myr non- 
accreting isochrone, as displayed in Figs, [I] and 
[3] As discussed in §3.1| calculations with a ^ 
cannot maintain luminous objects at ages of ~ 
Myr, with accretion rates consistent with the ob- 
served ones at this age, because of the too short 
thermal timescale. It thus seems difficult to ex- 
plain the luminosity of these objects with accre- 
tion. A possible explanation is that they are sig- 
nificantly younger (<§; 1 Myr) than the mean clus- 
ter age and experienced their episodic accretion 
phase quite recently. Another explanation, as dis- 
cussed in Chabrier et al. (2007), is that fast rota- 
tion and/or the presence of a magnetic field yields 
a smaller heat flux output, thus (i) a larger radius 
and (ii) a cooler T e g, while barely affecting the lu- 
minosity, for a given mass. The net effect would 
be to shift the location of ~ 1 Myr old objects 
of a given mass at cooler T e g for a given L, and 
thus on the right side of the 1 Myr non-accreting 
isochrone. Note that this may apply as well to 
the other (hotter) objects, even though there is no 
need to invoke such a process to reproduce their 
luminosity. This suggestion should motivate ob- 
servational determinations of the rotation velocity 
and the level of magnetic activity of these over- 
luminous objects. The present analysis supports 
the conclusion of Mohanty et al. (2009) on the 
origin of the T e g reversal in the young eclipsing 
binary brown dwarf 2M0535-05, which excludes 
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Fig. 3. — Evolutionary sequences in a HRD with 
episodic accretion. The labels close to each colored 
curve correspond to the cases described in Table 
[2] Symbols and other curves are the same as in 
Fig.0 



tion periods and projected radial velocities of low 
mass objects (Jeffries 2007). Episodic accretion 
thus seems to provide several matching pieces to 
the puzzle describing the early evolution of LMS 
and BDs. If our suggestion is correct, it has several 
drastic consequences: (1) what was interpreted as 
a significant age spread in SFRs or young clusters 
is essentially a spread in radius, thus in luminos- 
ity, for objects of comparable ages, (2) if young 
low-mass objects have experienced strong episodes 
of accretion during their embedded (class to I) 
phase, their contraction proceeds very differently 
from a standard, constant T e g, Hayashi track and 
still keep memories of these early episodes after 
about a few Myr, even if present accretion rates 
are negligible, (3) trying to infer the mass, thus 
an IMF, for young (< Myr old) clusters from 
mass-luminosity or mass-T c ff-Sp type relationships 
based on non-accreting objects/models very likely 
leads to results of low significance, and inferring 
the IMF for such young associations seems to be 
elusive. At this stage, determination of the proto- 
star/BD core mass function (CMF) with submil- 
limcter surveys will bring more robust information 
about the star formation process and the stellar 
mass spectrum (see e.g. Henncbellc & Chabrier 
2008, 2009). Finally, the present calculations show 
that the concept of a stellar birthline has no real 
significance, at least for low-mass objects, as the 
first appearance of these objects in a HR diagram 
is very random, due to the variety of prior accre- 
tion histories. 



an explanation based on prior accretion. At last, 
another explanation for the location of these ob- 
jects in the HR diagram is a less reliable photome- 
try. All the over-luminous objects in the Peterson 
et al. (2008) sample, for instance, are located in 
a region of the Orion Molecular Cloud with very 
high nebulosity (D. Peterson, priv. comm.). 

The suggestion, as explored in this Letter, that 
episodic accretion provides a plausible explanation 
for the observed spread in HR diagrams at ages of 
a few Myr (see Q, is supported by recent obser- 
vations of protostars (Enoch et al. 2009; Evans et 
al. 2009). Furthermore, our predicted significant 
spread in radius at ages of ~ Myr (see Fig. [2]) 
is consistent with the recently suggested existence 
of such a spread in the ONC, based on the rota- 
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